Sum Rule for the Optical Absorption of an Interacting 

Many-Polaron Gas 



J. Tempere and J. T. Devreese 
Departement Natuurkunde, Universiteit Antwerpen - UIA, Universiteitsplein 1, B2610 

Antwerpen, Belgium. 
(1/12/2000) 

Abstract 

A sum rule for the first frequency moment of the optical absorption of a 
many-polaron system is derived, taking into account many-body effects in 
the system of constituent charge carriers of the many-polaron system. In our 
expression for the sum rule, the electron-phonon coupling and the many-body 
effects in the electron (or hole) system formally decouple, so that the many- 
body effects can be treated to the desired level of approximation by the choice 
of the dynamical structure factor of the electron (hole) gas. We calculate 
correction factors to take into account both low and high experimental cutoff 
frequencies. 
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I. INTRODUCTION 



Sum rules have proven to be a powerful tool in the analysis of optical spectra |T|]. They 



provide expressions relating a frequency moment of the optical absorption spectrum to 
characteristic properties of a material such as its plasma frequency. Sum rules provide 
useful checks (on the optical properties) both for theory and for experiment. 

The goal of this paper is to derive a sum rule for the normalized first frequency moment 
of the optical absorption of a gas of continuum (Frohlich) polarons, including many-body 
effects between the charge carriers in this many-polaron gas. The normalized first frequency 
moment of the optical absorption is obtained from the optical conductivity Re[(7(co')] by 



This quantity is experimentally accessible: for example Calvani and co-workers [0 have 
determined (uj) as a function of doping (density) for the optical absorption bands attributed 
to polarons in a family of high-temperature cuprate superconductors 0. As such, the sum 
rule (|l]) for (u) applied to the polaron gas including many-body effects (such as interaction, 
screening, occupational effects of Fermi statistics,...) will prove a useful tool for the analysis 
of infrared spectra of such high-Tc materials. The study of infrared spectra in the framework 
of polaron theory can complement other indications of the presence of Frohlich polarons 
and bipolarons 0] in high-Tc materials, thus providing a more solid ground for hypotheses 
involving polarons and bipolarons in the mechanism of high-temperature superconductivity 



In the derivation presented here (in Section II), the many-body effects in the system 
of charge carriers are completely contained in the dynamical structure factor of the charge 
carrier gas, and they are formally factored out from the charge carrier - phonon interactions, 
which allows for a great deal of flexibility in treating the many-polaron gas. 
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We demonstrate how to circumvent two particular difficulties which complicate the prac- 
tical use of many-polaron sum rules in a comparison with experiment. The ffist difficulty is 
related to the presence of a delta function contribution at zero frequency and T = in the 



a sharp zero- frequency feature can be overlooked by experiment, and we derive a correc- 
tion factor -which turns out to be related to the density dependent effective polaron mass- 
to compensate for the inability to detect the delta function contribution. This correction 
factor is the subject of Section III. A second difficulty is that the theoretical many-polaron 
optical absorption spectrum contains a slowly decaying high-frequency tail which carries 
a substantial fraction of the spectral weight of the polaron optical absorption. In Section 
IV we discuss the influence of a high-frequency cutoff on the first frequency moment and 
propose a formula useful to extrapolate experimental data of the polaron optical absorption 
band at higher frequencies. 

II. FIRST FREQUENCY MOMENT OF THE OPTICAL ABSORPTION OF A 

MANY-POLARON GAS 

A. Exact expression 

The optical absorption coefficient is proportional to the real part of the optical conduc- 
tivity, Re[(T(c(j)]. To evaluate expression (|l]), both the zeroth and the ffist frequency moment 
of the real part of the optical conductivity have to be determined. The zeroth moment is 
easily related to the plasma frequency through the f-sum rule [H : 



In this expression, N is the number of polarons per unit volume, rrib is the band mass of 
the constituent electron (or hole) of the polaron, and e is the electron charge. To obtain the 
ffist frequency moment of Re[cr(co')], the optical conductivity is expressed as a current-current 



theoretical optical absorption spectrum of the polaron gas, at temperature zero 0. Such 
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correlation function according to the Kubo formalism of linear response theory (see e.g. 0). 
Integrating twice by parts, the current-current correlation function can be transformed into 
a force-force correlation function l|8|. This leads to: 



(3) 



where V is the volume (chosen equal to unity in what follows) and is the x-component 



of the force operator F = (i/h) -ff, X^^Li 



. To formulate the sum rule, the force- force 



correlation function is rewritten in spectral representation. For this purpose, we introduce a 
set {^^n} of basis states with corresponding energies {En}, of which \l/o is the ground state 
with ground state energy Eq. The spectral representation then becomes: 



e'-'{[F,{t),FM])dt = Y. e 



dt (4) 
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Denoting [En — EQ)/h as Uno, one arrives at the following expression for the first moment 
of the optical conductivity: 
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The result for the normalized first frequency moment of the optical conductivity is found 
from the expressions (|^) and (H): 
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To obtain a useful expression, the resummation over the unspecified set of basis functions 
should still be performed. This is done by using the integral representation 



— = lim III e'^^'^^'-^'dt 



(9) 



such that 
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where the Hermiticity of the force operator was used. If the response is time translational 
invariant (meaning that the response is only dependent on the time lapse since perturbation), 
the final result can be written as: 



lim 

Nhnih 6^+0 



(^o|i^.(t)i^.(0)|^o) e-''dt 



(12) 



Expression (|T^) is an exact expression, within the assumptions of linear response theory, 
for the sum rule for the normalized first frequency moment of the optical absorption of the 
many-polaron gas. 



B. The interacting Frohlich many-polaron gas 

Throughout this paper, we consider polarons in the continuum or Frohlich approximation 
0]. The system of interacting Frohlich polarons is then described by the Hamiltonian 

N 2 N N 2 / 

^ = E ^ + E ^^obtb, + Y: {VubuP. + VChipi) (13) 

]=1 k k ]=1 1^3=1 

where rj,pj are the position and momentum operators of charge carrier j with band mass 
TJib and change ibe, and aj^,ak are the creation and annihilation operators of a longitudinal 
optical (LO) phonon with frequency culo, Pk = expjik.rj} is the density operator 

of the charge carriers, Vk is the interaction amplitude matrix element characterizing the 
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interaction between the charge carriers and the LO phonons, and ^oo is the permittivity of 
the medium. The total force operator for the many-polaron gas, with both charge carrier - 
phonon interaction and Coulomb interaction between the charge carriers taken into account, 
takes the form 



I 

h 



N 



J2 k(^k&kPk - K^k Pk ; 



so that in the force-force correlation Frc{t)Fx{0), a factor |VkP oc a is present and formally 
factors out of the expression for {u) and Re[(T(co')] in equations (|12D and @, respectively. 
Herein lies an advantage of working with the force-force correlation as compared to the 
current-current correlation, especially at weak coupling. Thus, to calculate (ti;) (0) to lowest 
order in the charge carrier - phonon interaction strength a it is sufficient to evaluate the 
(a-independent) expectation value of the force-force correlation in a product wave function 
10) \{pe\) where \ipei) represents the ground-state wave function for charge carriers and |0) is 
the phonon vacuum: 
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The expectation value with respect to the phonon vacuum can be calculated to order a and 
results in 
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The remaining expectation value is the density-density Green's function G(k, t) 
-iQ{t) ((^ei |pk(t)Pk(0)| ^ci> /N so that 
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Introducing the dynamical structure factor S{k, u) as the spectral density function of the 
retarded density-density Green's function Gr of the electron (or hole) gas, we arrive at: 
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The modulus square of the Frohhch interaction amphtude is given by 
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where V (A) is the volume (surface). This leads to 
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The 2D and 3D results are related by the scaling relation 

(u;)2D («, 5*20) = (t^)3D (37ra/4, S'sd) (19) 

provided the corresponding 2D or 3D dynamical structure factor of the electron (hole) system 
is used. Taking the low-density limit of (0), the correct one-polaron scaling relation 



is retrieved. The present analysis shows that the generalization (|T^ of the one-polaron 
scaling relation holds in the many-polaron case, provided that the corresponding 2D or 3D 
dynamical structure factor is used. Both expressions (|19a|) and (|19b|) contain the Frohlich 
electron-phonon (or hole-phonon) coupling constant a only as a prefactor: the remaining 
integrations of ( |19a| ) and (|19b|) depend only on the details of the electron-electron (or hole- 
hole) many-body effects. 



C. Results and discussion 



For one polaron, the normalized first frequency moments are {uj)^^) ~ (2/3) (t^)2D 



(0^)20 = {Tc/2)aui^o [mi- These limits correspond to the low density limit of the expressions 



( p.9a|) and (|19b| ), which is taken by letting kp 0. In this hmit the dynamical structure 
factor ^(k, to) becomes strongly peaked along w = hk'^/{2miy). In the case of many po- 
larons, many-body effects such as the electron-electron interaction, screening effects, and 
the occupational effect due to Fermi statistics will play a role. 

Figure 1 shows the result obtained from ( |19aD ,( [W5D for the first frequency moment of the 
optical absorption in 3D and 2D as a function of density, expressed through the Fermi wave 
vector kp. For the 3D case, the material parameters of ZnO were used (culo = 73.27 meV, 
£00 = 4.0, TUh = 0.24 rrie) and for the 2D case, the material parameters for GaAs were taken 



(culo = 36.77 meV, £00 = 10.9, mb = 0.0657 rrie) [0. The limit of low density corresponds 
to the one polaron-results indicated by arrows starting on the y-axis. 

The dashed curve (for 2D) and the dotted curve (for 3D) were calculated with the 
Hartree-Fock structure factor of the electron (hole) system. The increase of (u) in the region 



fcp ~ 1 mfeC^LO / ^ shows an initial shift in spectral density towards higher frequencies. The 
full curve (for 2D) and the dash-dotted curve (for 3D) are the result obtained by substituting 
the RPA structure factor of the electron (hole) system in expressions ( |19aD ,( P!?5| ). When 
the many-body effects are taken into account on the level of the RPA, they result in a 
monotonously decreasing value of the normalized first frequency moment as a function of 
density. 

The results obtained from the sum rule expressions (|19a|) , (|19b|) are in agreement with 



the result ||T3| obtained by direct integration of the many-polaron optical absorption spec- 
trum calculated with the variational ground state wave function of Lemmens, Brosens and 



Devreese |jT^. Note that the derivation presented in the present paper does not rely on 
the variational wave function used by LDB. The use of the force-force correlation function 
introduces already a factor a in expression (p!^, so that it is sufficient to calculate the 
expectation value in (|1^ with respect to the phonon vacuum in order to obtain the first 
frequency moment up to order a in the charge carrier - phonon coupling strength. 



8 



III. ZEROTH FREQUENCY MOMENT AND MANY-BODY EFFECTS IN THE 

POLARON MASS 



In the polaron optical absorption spectra at zero temperature, a delta-function contribu- 
tion is present at = This feature of the polaron spectrum will, at low temperatures, 
elude experimental detection since all infrared absorption experiments have naturally a low- 
est measurable frequency. Nevertheless, this delta-function contribution will play a role in 
the zeroth frequency moment of the optical absorption and hence in all normalized frequency 
moments. Note that as the temperature is raised, this DC delta-function contribution will 



substantially widen and become detectable |]T5[. At temperature zero, we must rely on an 



adaptation of the f-sum rule for one polaron, derived in 

oo 

Nixe^ f ^ Nne 



+ / Re[a{u)]duj = . (20) 



2m* J 2mb 



In this expression, m* is the polaron effective mass. At very low temperature, the integral 
in the left-hand-side of (^) can be determined experimentally or used to derive the effective 



polaron mass from many-polaron optical absorption theories such as refs. |]T3|,|T6[. The LO 
phonon frequency cjlo is typically a few hundred cm~^, obviously within the range of far 
infrared spectroscopy. Expression (pOf) allows then to determine the polaron mass from the 
optical conductivity: 



oo 



m* Nire^ 



j Re[a{uo)]duj. (21) 



The evolution of the polaron effective mass as a function of the density of the interacting 
many-polaron gas is shown in Figure 2. To obtain the polaron effective mass, the sum rule 



was applied to the many-polaron optical absorption spectrum calculated in |13[ : 

2 



^^[^^^)] = i^^E^'l^'^l'^(k,^-.;Lo) (22) 



The result for {m* /mi, — 1)/q; is shown in Figure 2 as a function of density expressed via k-p- 
For low density the effective mass naturally converges to its one-polaron value in the weak 
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and intermediate electron-phonon coupling limit, given by m* = mb{l + a/6) in 3D and by 
m* = mb(l + 7ra/8) in 2D [|rT|. Both dashed curves were calculated using the Hartree-Fock 
structure factor. The shift of spectral weight towards higher frequencies, due to occupation 
effects, is reflected here as an initial increase in the effective polaron mass. The full and 
dash-dotted curves show the result using the RPA structure factor. The resulting effective 
polaron mass is a monotonously decreasing function of density. 

The normalized first frequency moment of the optical absorption, excluding the delta- 
function contribution, {uj)q, is accessible to experiment, using: 



u Re[a{uj)]du 

{u^)o~ • (23) 

J Re[a{uj)]duj 

The question of the high-frequency cutoff in the experiment will be discussed as a further 
refinement in the next section. The numerator in expression p3| ) is still coincident with the 
numerator in (|l]) for many-polaron optical conductivity. Based on (|2TD for the denominator, 
we obtain 

M=(i-^)Mo- (24) 

As a consequence, to obtain the normalized first frequency moment (uj) of the optical absorp- 
tion including the delta function contribution, one has to correct the result {uj)q obtained 

without the delta function, by multiplying it with a correction factor ( 1 ) . 

V m* / 

IV. HIGH FREQUENCY TAIL OF THE MANY-POLARON OPTICAL 

ABSORPTION SPECTRUM 

An upper bound uj^^^x of the integration domain (cutoff frequency) is necessarily present 
in the experimental determination of the frequency integrals in (0). Thus the experimen- 
tally accessible quantity for the experiment is 
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Re[a{uj)]du! 



(25) 



where Wmax is typically of the order of 10000 cm^^ 0. This upper bound will influence the 
experimental value of the first moment (uj) since a high-frequency tail {uj > c^max) is present 
in the optical absorption spectrum of polarons. 

Figure 3 illustrates this cutoff problem. The full curve in Fig. 3 shows the sum rule 
result for the normalized first frequency moment of the optical absorption of the many- 
polaron system as a function of the Fermi wave vector. The other curves show the results, 
obtained by direct integration of the spectrum calculated with the variational ground state 



wave function of ref. ||T^, using different values for the cutoff parameter cUmax- Already at 
i^max = 10 cjlo a considerable difference is found between the result with cutoff and the 
result without cutoff. 

For energy transfers hw much larger than the Hartree-Fock exchange energy of the many- 
polaron gas, the optical absorption will not be strongly modified by many-body effects. This 
property is used in the optical absorption measurements to extrapolate the high-frequency 
measurements of the reflectance by using a free-electron response (as described e.g. in ref. 
p!7|). For the many-polaron optical absorption in the high frequency limit (for example in 
3D) we find |13[: 



hm Ke[(T(co')J = -a- 



rrib 3 



(26) 



so that 



lim 

'-'max ^OO 



UJ iielaiuj)] = a 

rrib 3 



+ arcsin 



(27) 



For the first frequency moment this leads to 



lim 



(^)cxp (^max) + —a 



-|- arcsin 



(28) 
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Expression (pH]) allows one to correct the error made by using a cutoff frequency in determin- 
ing the experimental value for the first frequency moment of the polaron optical absorption. 
How large does uj^aax have to be for the limit ( pHf ) to be accurate ? Naturally, this depends 
on the density. The dynamical structure factor of the electron (hole) system reduces to the 
dynamical structure factor of a one-particle system for k ^ 2A;f, as does the static structure 
factor ||18|. At these values of the wave vector, the only frequency for which the dynamical 
structure factor S{k, uj — cjlo) differs substantially from zero is for uj — cjlo = A;^/2. Hence, 
we can estimate that the many-polaron optical absorption will be well approximated by the 
one-polaron optical absorption only for frequencies lo ^ 1k\^LO\,o. This is also the value of 
cUjnax which should be used as high-frequency cutoff in the experiments determining {oS)^^^. 
For example, in ZnO (culo = 73.27 meV, = 4.0, = 0.24 mg) with a charge carrier 
density of n = 10^° cm~^, this would correspond to cjmax ~ 18.9 cjlo or ^^max ~ 11000 cm~^. 

V. CONCLUSION 

Expressing the optical conductivity as a force-force correlation function, a closed expres- 
sion is obtained for the sum rule of the normalized first moment of a many-polaron system, 
in such a way that the electron-phonon coupling and the many-body effects of the electron 
(hole) system formally decouple. This procedure allows one to treat the many-body effects 
in the system of charge carriers using any desired approximation (Hartree-Fock, RPA, ...). 
The results obtained by the sum rule derived here are valid for a many-polaron system with 
weak electron-phonon coupling and at temperature zero. 

Experimental data are characterized by a cutoff frequency both for low and high frequen- 
cies This in general complicates the use of any sum rule, since its evaluation in principle 
requires a summation over all frequencies, lo = O...00. To overcome this problem we in- 
troduced formula ( pSf ) for many-polaron optical absorption spectra at high frequencies, and 
discussed up to which frequency experimental data have to be available so that this extrap- 
olation formula is useful. At zero temperature, a delta function contribution is present at 
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c<j = in the many-polaron optical absorption spectrum 0. The spectral weight of this delta 
function contribution is related to the polaron mass which in its turn can be derived from 
a measurement of the zeroth frequency moment of the many-polaron optical absorption. 
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FIGURE CAPTIONS 

Figure 1: The normalized first frequency moment of the optical absorption of a gas 
of polarons including many-body effects in the system of constituent charge carriers, as 
obtained with the present treatment, is shown as a function of density, expressed through 
the Fermi wave vector fcp in polaron units. The top two curves are for the 2D polaron gas 
(with parameters corresponding to GaAs: cjlo = 36.77 meV, e^o = 10.9, m?, = 0.0657 m^), 
and the bottom two curves are the results for the 3D case (with parameters representative 
of ZnO: ujlo = 73.27 meV, = 4.0, rub = 0.24 rUe). 

Figure 2: The effective polaron mass in the many-polaron system is shown as a function 
of density, for the 2D case (GaAs parameters: culo = 36.77 meV, ^oo = 10.9, nib = 0.0657 
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me) and the 3D case (ZnO parameters: culo = 73.27 meV, £00 = 4.0, rrib = 0.24 me), and 
for different approximations (Hartree-Fock and RPA) to take tlie many-body effects into 
account. 

Figure 3: The effect of introducing an upper hmit (a cutoff, cumax) to the integration 
domain in the expression for the first frequency moment (u) of the optical absorption of the 
polaron gas is shown in this figure. The different curves show {u) / {uj)f^^^Q as a function of 
the Fermi wave vector for the set of cutoff values listed in the legend, obtained from direct 
integration of the many-polaron spectrum as calculated in [|^. The full curve shows the 
sum rule result (without cutoff) in the present treatment. In the inset, the absolute value 
of (uj) is shown as a function of the Fermi wave vector, for the same set of cutoff values. 
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2D polaron gas, Hartree-Fock 



-^7C/2 



-^2/3 



2D polaron gas, RPA 
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